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1 Introduction 


The computation of the expectation (r^) of the electron for atoms with a single electron is 
a standard problem in quantum mechanics, see [19, 29]. Eor a given energy state n, the problem 
is expressed as 

COO 

Jo 

where Rnt{r) is the radial solution of the Schrodinger equation for the hydrogen atom. Condi¬ 
tions on the parameters n, I, k are determined by the convergence of this integral. 

In the non-relativistic situation, the solution is given in terms of the Hahn polynomials [4]: 


ht^\x,N) 


(1 - iV)m(/3 + l)m [ -m, a +(5+ m + l, -X 

d + l.l-N 



( 1 . 1 ) 
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In particular, these expectations are given in terms of the Chebyshev polynomials of discrete 
variables [18, 24] 

tm{x,N) = AI) 


in the form 

= ^(2/u)“^4+i(n-£-1,-2^-1), when /c =-1, 0,1, 2,..., (1.2) 

and 

(r^) = —(2^)“^t_fc_2(n-£-l,-2£-1), when k = -2,-3,...,-21 - 2. (1.3) 

^ ' 2n 

The parameters are /x = Z/na^ with oq = f? jra^ the Bohr radius and Z is the nuclear charge. 
The constants m and e are the mass and charge of the electron, respectively. 

The identity 


^ _ 1 , _ 2 ^ _ 1 ) = 

, ) r(2£ + 2) 

follows from (1.1). Then (1.2) becomes 


—k, k + 1, —n + £ + 1 
l,2i + 2 


/ k\ 1 r(2t' + /c + 3) ^ k — 1 — k, k 2, —n 1 

y /n£- 2n{2n)’^ r{2i + 2) ^ \ 1,21 + 2 


for k = —1,0,1,2,..., and (1.3) 




1 T{2l-k) 

2n{2^if T{2i + 2) 


3-^2 


—2 + k, —k + 3, —n + 1+1 

1,21 + 2 


for k = —2, —3,..., —21' — 2, where the dependence upon the parameters n and ^ have been 
made explicit. 

In the quantum case, the radial component of the wave function for a hydrogen atom with 
nuclear charge Z is characterized by two quantum numbers: n the principal quantum number 
and i the orbital number. The corresponding normalized radial function is 


Rni{r) = Ani{2+rf eyip{-pLr)L^^^l_^{2pLr), 


where the normalization constant is 


An£ 


(2/i)3 (n — I — 1)! 
2n {n + £)l 


and 


(x) = 


r(a + m + 1) 


r(m + l)r(l + a)^ ^ 1,^1 + q; 


—m 


is the associated Laguerre function; see [14, formula 8.972.1]. The expectation value of a power 
of the radial coordinate is given by 


roo 

Jo 

with neN, 0<l<n — 1, and k £ 7,. 


(1.4) 
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The goal of the work is to compute the integral in (1.4) by the method of brackets, to illustrate 
its flexibility. The reader will find in [3, 6, 9, 10, 11, 12, 16] a collection of examples of definite 
integrals evaluated by this method. The basic procedure is described in Section 3. 

The examples presented here are to be considered as the beginning of a series of calculations 
of integrals related to the hydrogen atom. These include the evaluation of the integral [13] 

poo 

JnL = / 

Jo 

given by S.K. Suslov and B. Trey [25]. The method of brackets provides an alternative method 
of proof that only uses the hypergeometric representation of the Laguerre function. The method 
can also be used to discuss the relativistic situation. Details will appear elsewhere. 

The reductions of the formulas discussed here uses basic properties of the gamma function, 
such as 


r(a + n) = r(a)(a)n and 


(o)—n 


(-ir 

(1 - a)n 


for a G M, n G N. 


(1.5) 


Here {a)n = a{a + 1) • • • (a + n — 1) is the Pochhammer symbol. 


2 A direct evaluation 


This section presents a direct evaluation of the integral 

poo 

= {2^rAl, / (2.1) 

Jo 

given in (1.4). The proof is based on some identities for the associated Laguerre function 
appearing in the integrand. The methods presented here are then compared with the evaluation 
by the method of brackets explained in the next section. 

The first identity used to modify the integrand appears in [14, formula 8.976.3] 


Therefore 


2 _ r(Q! + m + 1) ^ f 2m — 2s\ r(2s + 1) 


■E 


22”^r(m + 1) ^ y m — s J r(a + s + l)r(s + 1) 


Lifi2x). 


(r^\ = i2uS^^A^ r(£ + n + 1) 

\' /ni ^’^^22(^-^-i)r(n — 7) 


n—£—1 

E 

s=0 


2 (n - 7 - 1 - s)\ r( 2 s + l) 

n - £ - 1 - s y r(2£ + 2 + s)T{s + 1) 


where 


poo 


To obtain an expression for Gi^k,s{h'), the representation 


—n 


ra.( \ ^ r(a + n + 1) 

r(n + i)r(i + a)^ ^Vi + 


X 


( 2 . 2 ) 


(2.3) 


(2.4) 


(2.5) 


for the Laguerre function (see [14, formula 8.972.1]) is used. 
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Theorem 2.1. The integral Gi^k,s{h) given by 


^ . , r(4£ + 2s + 3)r(2^ + A: + 3) ^ 

Ge,k,s{h) = 7^77;— , 2J1 


r( 2 s + l)r( 4 £ + 3 )( 2 /i) 2^+^+3 

Proof. The hypergeometric representation (2.5) shows that 


-2s, 2f+ A: + 3 
4^ + 3 


. 2 ( 2 £+i)., ^ r(4f + 2s + 3) / 

^ ^ ^ r(2s + i)r(4f + 3) ^ ^ V 


-2s 
4f+ 3 


4/ir 


Expanding the hypergeometric function gives 
r(4^ + 2s + 3) 


G£^k,s{h) — 


/■“ ^ (-2s)j {A^iry oe+2+k -2u.r, 

r( 2 s + l)r( 4 A + 3 )yo ,t^( 4 ^ + 3 ), j! 


r(4f + 2s + 3) ^ • /■“ 2 i+ 2 +k+j - 2 £Lr^ 

r(2s + l)r( 4 £ + 3) {U + 3)jj\^ Jo 


r(4f + 2s + 3) 


2s 

E 


(-2s)j(4/r)Jr(2^ + A: + j + 3) 


r(2s + l)r(4A + 3) ^ (4f + 3)jj! (2/x)2^+fc+i+3 


r(4f + 2s + 3) 


2s 

E 


{-2s)j2i 


r(2s + l)r(4£ + 3)(2/i)2^+fc+3 ^ (4f + 3)jjl 


—r(2f + /c + j’ + 3) 


r(4f + 2s + 3)r(2A + A: + 3) (-2s)j(2f + A: + 3)j 


r(2s + l)r(4A + 3)(2;u)2<?+fc+3 


E 

j=0 


{M + 3),j\ 


-2^' 


r(4f + 2s + 3)r(2A + k + 3) ^ (-2s)j(2f + A: + 3)j 


r(2s + l)r(4£ + 3)(2/i)2^+fc+3 

r(4f + 2s + 3)r(2A + /s + 3) 
“ r(2s + l)r(4A + 3)(2^)2«+fc+3 

This is the stated form for G£^k,s{h)- 


E- 

j=0 

2^1 


iM + 3)jj\ 

-2s, 2A + A: + 3 
4^ + 3 


-2^ 


Note 2.2. Observe that s G N, so the hypergeometric function in the expression for G£^k,s{h) 
is actually a polynomial in its last variable. Thus, there are no convergence issues. 

The expression for Gi^k,s{h) (2.3) are used to produce the next result (after the change 
SI—)>n — f—1 — s). 


Corollary 2.3. For n = 1,2,..., i = 0,1,..., n — 1 and A: G Z with 2i + k + 3 > 0. The 
moments of the hydrogen atom are given by 



T{2l + k + 3){2n + 2l)\ 
n22r^-2^-i(4£ + 2)!(2/i)*^(n + f)!(n - f - 1)! 


n-^-l /n+t\(n-£-l\ 

\ ^ V s / V s ) TP 

^ — J2FTM\ —2hi 

s=0 V 2s / 



f-l-s), 2A + A: + 3 
4^ + 3 



Note 2.4. The restriction 2A'+A;+3 > 0 avoids the singularities of the gamma factor r(2A'4-A;+3). 
Also observe that the first entry in the series 2 E 1 in the answer is a negative integer, therefore 
the series reduces to a finite sum. 
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In this article the expectation values of the powers of the radial coordinate of the hydrogen 
atom in a framework of quantum mechanics, that is, in the non-relativistic case are computed. In 
the Introduction it was stated that this already has appeared in the literature. In the relativistic 
case, results for these expectation values of the powers of the radial coordinate appeared in 2009. 
Indeed, the relativistic Coulomb integrals are contained in [22, 23]. The treatment of the results 
obtained in [22, 23] by computer algebra methods is described in [17, 20]. 

In the nonrelativistic case of quantum mechanics, the corresponding questions were success¬ 
fully solved by direct calculation. For example, in [21] useful relations between different Laguerre 
polynomials were found. In [5] the radial expectation values are given for D-dimensional hydro- 
genic states with D > 1. The same quantities are discussed in a more general setting in [26, 27]. 
The radial expectation values of hydrogenic states in momentum space appear in [28], repre¬ 
sented in terms of Gegenbauer polynomials instead of Laguerre polynomials. All these results 
were obtained by direct calculations too. The method of brackets may significantly simplify the 
calculations for these tasks. This will be discussed in a future publication. 

The method of brackets is not the unique successful method which involves integral transfor¬ 
mations. Traditional methods based on Mellin-Barnes transformation may be efficient tools in 
order to obtain new results in quantum field theory [1, 2, 7, 8, 15]. 


3 The method of brackets 


The evaluation of the integral giving the mean value (r^) (2.1) presented in the previous section, 
used the relation (2.2) in a fundamental way. A method to evaluate integrals over the half line 
[0, oo), based on a small number of rules has been developed in [11, 12]. This method of brackets 
is described next. The heuristic rules are currently being placed on solid ground [3]. The reader 
will find in [6, 9, 10] a large collection of evaluations of definite integrals that illustrate the 
power and flexibility of this method. 

For a G C, the symbol 

roo 

{a) = / x°'~^dx 

Jo 

is the bracket associated to the (divergent) integral on the right. The symbol 


'fn 


(-ir 

F(n -I- 1) 


is called the indicator associated to the index n. The notation or simply denotes 

the product ''' 4>ir- 

Rules for the production of bracket series. 

Rule Pi. Assign to the integral f{x) dx a bracket series: 


4>na{n){an + fd). 

n 


Here the coefficients a(n) come from an assumed expansion f{x) = ^ c/)na{n)x'^"'^^ The 

n>0 

extra ‘—1’ in the exponent is set for convenience. The coefficients are written as a(n) because 
these will soon be evaluated at complex numbers n, not necessarily positive integers. Now we 
need to state how to convert the bracket series into a number. 

Rule P 2 . For a G C, the multinomial power (ai-|-a 2 -|-- • -- 1 - 0 ,.)“ is assigned the r-dimensional 
bracket series 


ni 712 yir 


nr -t- ni H-h 

F(-a) 
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Rules for the evaluation of a bracket series. 

Rule El. The one-dimensional bracket series is assigned the value 

'^4>nf{n){an + b) = -^f{n*)T{-n*), 

n ' ' 

where n* is obtained from the vanishing of the bracket; that is, n* solves an + b = 0. This is 
precisely the Ramanujan’s master theorem. 

The next rule provides a value for multi-dimensional bracket series of index 0, that is, the 
number of sums is equal to the number of brackets. 

Rule E 2 . Assume the matrix A = (atj) is non-singular, then the assignment is 

E-E (pni-'-Ur /(ni, . . . , nr-)(aiini H-h airUr -|- Cl) • • • (Or-ini H-h OrrUr -I- Cr) 

ni Ur 

= • • • ’ • • • r(-<), 

where {n*} is the (unique) solution of the linear system obtained from the vanishing of the 
brackets. There is no assignment if A is singular. 

Rule E 3 . Each representation of an integral by a bracket series has associated an index of 
the representation via 

index = number of sums — number of brackets. 

It is important to observe that the index is attached to a specific representation of the integral 
and not just to integral itself. The experience obtained by the authors using this method suggests 
that, among all representations of an integral as a bracket series, the one with minimal index 
should be chosen. 

The value of a multi-dimensional bracket series of positive index is obtained by computing 
all the contributions of maximal rank by Rule E 2 . These contributions to the integral appear 
as series in the free parameters. Series converging in a common region are added and divergent 
series are discarded. Any series producing a non-real contribution is also discarded. There is no 
assignment to a bracket series of negative index. 


4 The evaluation of the expectations. A first bracket calculation 


This section describes the evaluation of the integral 

/*oo 

In,e,k{h)--= / ^(2H] dr, 

Jo 


(4.1) 


that appeared in (1.4) by the method of brackets. The expectation value of a power of the radial 
coordinate is then given by 

This integral can be scaled to 
1 


Jn,£,kih) — 




(2/i)3+2^+^ Jq 

This does not appear in the table [14]. The closest entry is 7.414.10: 


(4.2) 


n [ L“(x)l ^dx = 2)r(^+ D 

Jo 7r(n!)262a+i 


r(a -I- n -I- 1 ) 2E1 


—n, a -b 5 


— n 


2\2 


(1-D 
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Note 4.1. In the evaluation of (4.1), it is convenient to write it as 


I, 


poo 

{^^) := / r2+2^+^e-^^L2i+i_i(Hr)Lf_+/_,(Cr)dr 

Jo 


n,e,k:A,B,C\IJ‘) ■ = 

and then consider the limiting value as A, B, C tend to 2fj.. 

The computation of (4.2) described in this section is obtained without any further identities 
for the Laguerre function. Next section describes the computation of the function Gi^k,s{p)i 
defined in (2.4). 

The first step is to compute a series representation for the factors in the integrand. 

Lemma 4.2. The functions in the integrand of (4.1) have series given by 

e-“" = 


ni 


and 


LZiix) = r(a + 1 + m) ^ 

712 ^ 


X 


.712 


+ m- n 2 )r(l + a + n 2 )' 


Proof. The series of the exponential function is elementary. Indeed, 


g-ax ^ 


E 

ni>0 


(-“) 


ni 


= 


ni! 


E 

ni>0 


(- 1 ) 


ni 


TT-l! 






To evaluate the series of the Laguerre function, treat m as a real non-integer parameter, and 
observe that 

r(a-M-hm) ^ {-m)n2 

■^m\X) -n/^. I lA-n/.™ i i\ / . 


r(a -I- l)r(m -I- 1) (a + 1)^2 ^ 2 ! 


r(a -|- 1 -I- m) 


E 


r(n 2 — m) 


,.n 2 


r(m-Fl) r(-m)r(a-h 1 + n2) 712! ’ 

The series for the Laguerre function now follows from the identity 


r(n 2 — m) 


= (-l) 


712 


r(i -|- m) 


r(—m) r(l-|-m —77,2) 

valid for 772 G N and ttt, 0 N. ■ 

The series given in Lemma 4.2 are now used directly to evaluate the integral (4.1). This gives 


T 




r 

{h) = 

Jo 


,2+2l+k 




L Til 


E 


L ^2 


E 


L ^3 


r(£ + n + l) 


r (77 - ^ - ?T,2)r(2f + 2 + n 2 ) 




n2 


r(£ + n + l) 


r (77 - ^ - ?T, 3 )r( 2 f + 2 + n^) 




713 


dr 


poo 

/ ^2+27+fc+ni+n2+n3^^^ni^n2(^n3^^^^ 

Jo 


712,713 


711,712,713 
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r 2 (£ + n + 1 ) 

^ r(n - n2)r(2^ + 2 + n2)r(n - n3)r(2^ + 2 + ns) 

= {ni + n2 + U'i + ^ + 2^ + (l)ni,n2,ni 

ni,n2,nz 

r 2 (£ + n + 1 ) 

^ r(n - n2)T{2(. + 2 + n2)r(n - I- n^)T{ 2 l + 2 + ns)' 

This intermediate result is stated next. 


Proposition 4.3. A bracket series for the integral In,£,k-,A,B,c{h) given by 
In,i,k;A,B,c{h) — (^1 + ^-2 + ns + 3 + 2.^ + k) A^^ <f)rii,n 2 ,n^ 

ni,n2,nz 

r2(£ + n + 1) 

r(n — a. — n 2 )T{ 2 ^ + 2 + n 2 )r(n — i — ns)r(2£ + 2 + ns) 

The bracket series above contains one bracket and three indices, thus it is expected that the 
method will produce a double series as an expression for In,e,k-,A,B,cih)- The vanishing of the 
bracket gives 

ni + n2 + ns = —3 — 2i — k, (4.3) 


showing the two free indices. 

Solving for ns. Replacing ns = —ni — n2 —t, with t = 2i + k + 3, in the bracket series yields 
the expression 


I, 


n,£,k-,A,B,C 


(h) = 


r 2 (n + £ + l) 


OO 


r(ns+n 2 + t) (-§) 


B\n2 


niT^=o ^ - n2)T{2i + 2 + n2)r(ni + n2 + s)r(-l - k - n\ - n2)ni\n2\ 


with s = n + ^ + 3 + A;. Using (1.5) yields 

^ r^(n + £ + i)r(t) 

n,i,k-,A,B,c{h) - £)r{2e + 2 )r(s)r(-i - k) 

{t)ni+n2 (1 - n + i)n2 (fc + 2)ni+n2 / . xn-;> f 

{2i + 2)n2{s)ni+n2ni\n2l \C J \C J 


Then use 


(^)ni+n2 — (5)n2 (^ T n2)ni 
to produce 

. . r^(n + £ + i ) r ( t ) 

n,e,k-,A,B,c{h) - £)r{ 2 £ + 2) r ( s ) r(-l - k) 

{t)n2 (^ + ^2)ni (1 — n + £)n2 {k + 2)^2 (fc + 2 + n2)ni / . .^^2 

(2£ + 2 )„ 2 (s)n 2 (s + f^ 2 )nini!n 2 ! \C J [cj 

The sum corresponding to the index m, which appears only in 3 places, is chosen as the internal 
sum. This yields 

r fn)- r^(n + £ + i)r(t) 

n,£,k;A,B,C[hJ + 2 )r(s)r(-l - k) 
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xE 

n 2=0 


(i)n2 {k + 2)^2 (1 — n + 


)n2 


B\ ^ (t + n 2 )ni(^ + 2 + 712) 

nT^O ^ 


ni 


(2£ +2)n2(s)n2«2! V Cj {s + n2)nini\ 

The inner sum is now identified as a hypergeometric function to produce 

r 2 (n + £ + i)r(t) 


ni 




xE 

n 2=0 


C'*r(n - i)r{ 2 i + 2 )r(s)r(-i - k) 

it)n 2 i (k + 2)^2 (1 — n -\- k)n 2 


{2i + 2)„2 ('S)n 2 ^ 2 ! 


/t + 712, 2 +/c + n 2 

-C s + m_ 


Note 4.4. The same procedure can be used to treat the cases obtained by solving for tii or 772 
in the equation (4.3). The corresponding integrals are 

T^n + £ + l)r{t) 


r 


(1) 


n,t,k-,A,B,C 


(^) = 


A*r 2 (n-£)r 2 ( 2 £ + 2 ) 

(t)n 2 (l ~ 71 + i)n2 


E 

n 2=0 


B 


(2£ +2)^2 712! Vx4 


712 


2 F 1 


t + 7i2, 1 — n + £ 
2 £ + 2 


and 


r 


( 2 ) 


n,e,k-,A,B,C 


(^) = 


r 2 (n + £ + i)r(t) 
A*r2(n-£)r2(2£ + 2) 


E 

n3=0 


( 2 £ + 2 )^ 3713 ! 


A 


713 


(t)n3,(l 71 + £)n3 


t + n 3 , 1 — n + £ 

2 £ + 2 


At this point, the parameters A, B, C are replaced by the value 2/i, in order to continue the 
evaluation. This gives 

r 2 (n + ^ + i)r(t) 


7n,£,k{lF) 


( 2 /i)*r(n - £)r{ 2 £ + 2 )r(s)r(-i - k) 

2 F 1 


xE 

^2=0 


it)n2 (k + 2)n2 (1 — 71 + £)n2 / ^ ^no tt^ f t + n2, 2 + k + n2 


{2£ + 2)jj 2 (s)n2''2'2! 


5 + 712 


Observe that 1 — n + f is a negative integer, so this is actually a finite sum. Using Gauss’ 
evaluation 


2 i"l 


a, b 


1 = 


r(c)r(c — a — b) 
r(c — o)r(c — b) 


for c — a — 6 > 0, 


and expressing the resulting gamma factors in terms of Pochhammer symbols to obtain 

r( 7 i + £ + l)r( 2 £ + k + 3)r(n -i-k-2) 


7n,£,k{k'} — 


( 2 / 7 ) 2 ^+fc+ 3 r 2 (n - £)r{ 2 £ + 2 )r(-l - k) 

{k + 2)n2{^ — n F £)n2i‘2£ F k + 3)r, 

n2=0 


(2£ + 2)^2(£ F k F 3 — n)n2'^2^- 


The final step identifies this series as a hypergeometric series to produce: 
r(7i + £ + l)r(2£ FkF 3)r(n -£-k-2) 


7n,£,k{l^) 


(2/x)2^+fc+3r2(n - £)r(2£ F 2)r(-l - k) 
k F2, 1 f£ — n, 2£FkF3 
2£F2, lFkF3-n 


X 3 F 2 


The results of this section are summarized in the next statement. 
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Theorem 4.5. For n, i, k as above, 



r{2i + k + 3)r{n-£-k-2) 

2n(2^)^r(n - £)T{ 2 i + 2)r(-l - k) ^ ^ 


kk + 2,l + i-n,2i + k + 3 
\ 2i + 2,l + k + 3 — n 



5 The evaluation of the expectations. A second approach 

The moment {r^)n(. has been expressed in (2.3) as a finite sum values of the integral 

COO 

+')(4^r)dr, (5.1) 

Jo 

where the index s is an integer varying from 0 to n — I’ — 1. Corollary 2.3 provides an expres¬ 
sion for {r^)ni as a finite sum of values of the hypergeometric function 2 F 1 evaluated at the 
argument 2. The hypergeometric terms appearing in the mentioned representation are actually 
finite sums, so the convergence of the series is not an issue. An alternative form is derived in 
this section that extends the range of validity of Gi^k,s{k) to a larger range for the parameter s. 

The goal is to produce a representation of the series for the Laguerre polynomials, given 
initially by 


l"iri = h(a -b n -H) / 

r(n + l)r(a + l)^ 


—n 
a ^ 


This series is now written in a form suitable for the application of the method of brackets: 

r(a-bn-M) ^ 


r(n -|- l)r(Q; -|- 1 ) + t)A;i ^1! 


r(Q; -|- n -|- 1 ) / 1 \ / \ 

r(n + l)r(a + l) (-^)fci(-a)-fci 


A:i=0 




r(Q; -|- n -|- 1 ) , u, 

r(,. + i)r(a + i) ^ 

ki 

p/ I nrr -Li-irr - Irv-V </>ir -n + fci r(-a - ki)x 

r(n -I- l)r(a -I- l)r(—n)r(—a) ^ 


fcl 


To produce a bracket series representation of the last expression, observe that 




i 


and this leads to 


Ln{x) 


r(a -|- n + 1) 

r(n -|- l)r(a -|- i)r(—n)r(—a) 


(t)i2z{-n + ki + k2){-a 


kl + k3)x^F 


The vanishing of the brackets provides two representations for the Laguerre function, denoted 
by Tj. 

Case 1. Take ki as a free index. Then k^ = n — ki and k^ = ki + a yields the expression 


r(Q;-|-n + l) ^ {-n)ki x^^ 
r(n + l)r(a 1 ) (a + 1 )^^ ki\ ' 


This is the original series for L“(x). 
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Case 2. Take k 2 as a free index. Then k\ = n — k 2 and k^ = a + n — k 2 yields 


T2 = 


r(a + n + l)x^ 


r(n + l)r(a + i)r(—n)r(—a) 


r(-n + k 2 )T{-a - n + /C 2 ) 


-X 


1—^2 


k2=0 


k2\ 


(5.2) 


Case 3. Taking k^ as a free index does not produce a representation for T"(x). 

The next step is to use the T 2 representation to evaluate the integral Gi^k,s{k)- By equa¬ 
tion (5.2), the expression for T“(x) is now written as 


L“(x) 


r(a + n + l)x” 

r(n + l)r(a + l)r(-n)r(-a) 


+j)T{-a-n +j)x C 

j=0 


Using this representation in (5.1) produces 

r(4£ + 3 + 2s)(4^)2" 




r(2s + i)r{M + 3)r(-2s)r(-4^ - 2 ) 

00 ^CX) 

V ^jr{-2s + j)r{-Ae - 2 - 2s + j)(4/i)“^' / 

Jo 


j=0 


Evaluating the last integral in terms of the gamma function and simplifying produces a proof of 
the next result. 


Theorem 5.1. The integral 

/*oo 

+')(4^r)dr 

Jo 

is given by 

4" r(3 + 2£ + A: + 2s) /-2s, -2s - 4£ - 2 1\ 

r(2s-hl) ^ -2-2i-k-2s 2)' 


6 A couple of examples 


The method of brackets has been used here to produce analytic expressions for the mean radius 

/•OO 

Jo 

stated first in (1.4). The physically relevant parameters are 


n = 0,1,2,..., 0 < / < n — 1, fc G M. 


The expressions include 

r(2/ +A: + 3)(2n + 2/)! 




^ 22 n- 2 £-i(' 4 £ + 2)!(2/i)*^(n -h /)!(n - / - 1)! 


n—l—l (n+ll\ /n—£—l 


E 

s=0 


/2n-|-2A 
f 2s / 


>Ei 


-2(n-/-l-s), 2e + k + 3 
U + 3 


( 6 . 1 ) 


where {r^)ni is given as a finite sum of hypergeometric terms and 




r(2/ + A: + 3)r(n-/-A:-2) 
2 n( 2 ^)fcr(n - /)r( 2 / + 2 )r(-l - k) 


3 F 2 


k + 2, 1 + i — n, 2£ + k + 3 
2£ + 2, l + k + 3 — n 


given in Theorem 4.5. This section compares these expressions with the results found in the 
literature. 
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Example 6.1. Take £ = n — 1. Then the sum (6.1) reduces to 1 since the index s must vanish. 


Then 



In particular, for A; G N, this becomes 



Example 6.2. The case £ = n — 2 reduces the sum (6.1) to two terms. The result is 



{k‘^ + 3k + 2n)r{k + 2re — 1) 
2(2^)^(2n- 2)! 
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